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GFOM£TP.Y AND STARVATION EFFECTS IN HYDRODYNAMIC LUBRICATION 

by David Brewe 
Propulsion Laboratory 

AVRADCOM Research and Technology Laboratories 
Lewis Resea-ch Center 
Cleveland, Ohio 

and 

Bernard J. Hamrock 
Lewis Research Center 
Cleveland, Ohio 


SUMMARY 


I 

u 


Numerical methods were used to determine the effects of lubricant starvation on the 
miinimum film thickness under conditions of a hyarody namic point contact. Starvation was 
effected by varying the fluid inlet level. The Reynolds boundary conditions were 
applied at the cavitation boundary and zero pr'^ssure was stipulated at the meniscus or 
inlet boundary. A minimum-f ilm-thickness equation as a function of both the ratio of 
dimensionless load to dimensionless speed and inlet supply level was determaneo. By 
comparing the film generated under the starved inlet condition with the film generated 
frorr the fully flooded inlet, an expression for the film reduction factor was obtaineo. 
Based on this factor a starvation threshold was defined as well as a critically starvec 
inlet. The changes in the inlet pressure buildup due to changing the available lubri- 
cant supply are presented in the form of three-dimensional isometric plots and also in 
the form of contour plots. 


NOMENCLATURE 


"o'^l 

D 

H 


in 

IT 

in 


H' 


h 

I. 

N 

P 

P 

R 

W/U 

u 

w 

X 

X 

Y 
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a 

B 


<p 

Subscripts 
cr ■ 

f 

x,y ■ 


least-squares coefficients 
difference | (H^ - H^I/Hq’xIOO, percent 
dimensionless film thickn'ss, h/R^^ 
dimensionless minimum (cent'al) film thickness, 
dimensionless calculated minimum (central) film thickness 
dimensionless fluid inlet level, h^^/R^ 
dimensionless fluid inlet level (onset of starvation) 
film thickness, mi 

minimum (central) film thickness, m 
reduced hydrodynamic lift, dimensionless 
direction normal to boundary 
dimensionless pressure, pR /n u 
2 

pressure, N/m 


* effective radius of curvature. 




+ R„) 


= ratio of dimiensionless load to dimensionless spee'.; 

•= average surface velocity in x direction, (u^ -*• Ug)/2, m/s 

* load capacity, N 

* dimensionless coordinate, x/R^^ 

« coordinate along rolling direction, mi 
» dimensionless coordinate, y/R^ 

» coordinate transverse to rolling direction, m 
■ radius ratio, Ry/Rx 

* film reduction factor 

* fluid viscosity at standard temperature and pressure, Ns/m^ 
Archard-Cowking side-leakage factor, 1/(1 + 2/(3 q)) 

critical 

flooded conjunction 
coordinate direction 
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GEOMtTPy AND STARVATION EFFECTS IN HYDRODYNAMIC LUBRICATION 

by David Brewe 
Propulsion Laboratory 

AVRADCOM Research and Technology Laboratories 
Lewis Research Center 
Cleveland, Ohio 

and 

Bernard J. Hamrock 
Lewis Research Center 
Cleveland, Ohio 


SUMMARY 


I 

u 


Numerical methods were used to determine the effects of lubricant starvation on the 
minimuni film thickness under conditions of a hyoroaynamic point contact. Starvation was 
effected by varying the fluid inlet level. The Reynolds boundary conditions were 
applied at the cavitation boundary and zero pressure was stipulated at the meniscus or 
inlet boundary. A minimum- f i Im- thickness equation as a function of both the ratio of 
dimensionless load to dimensionless speed and inlet supply level was determineo. By 
comparing the film generated under the starved inlet condition with the film generated 
from the fully flooded inlet, an expression for the film reduction factor was obtained. 
Based on this factor a starvation threshold was defined as well as a critically starvec 
inlet. The changes in the inlet pressure buildup due to changing the available lubri- 
cant supply are presented in the form of three-dimensional isometric plots and also in 
the form of contour plots. 
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pressure, pP 

/n u 
x' 0 


pressure, N/m 

effective radius of curvature. 


F R / (R + 
X y^ ' X 




= ratio of dimensionless load to dimensionless speeo 


Ug)/2, m/£ 


average surface velocity in x direction, (u^ + 
load capacity, N 
dimensionless coordinate, 
coordinate along rolling direction, mi 
dimensionless coordinate, y/Rj^ 
coordinate transverse to rolling direction, m 
radius ratio, Ry/Ry 
film reduction factor 

fluid viscosity at standard temperature and pressure, 
Archard-Cowking side-leakage factor, 1/(1 + 2/(3a)) 

critical 

flooded conjunction 
coordinate direction 


Ns/m 
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18-2 INTBODl'CTION OF POOR QUALITY 

The effect of starvation in a hydrodynamical ly lubricated conjunction can be 
studied by systematically reducing the inlet supply and observing the resultant pressure 
distribution and film thickness. This starvation effect can have a significant role in 
the operation of machine elements. For example, roller-end wear due to roller skewing 
can be a critical problem for high-speed cylindrical roller bearings. It is desirable 
that the hydrodynamic film generated between the roller end and the guide flange provioe 
stiffness and damping to limit the amplitude of the roller skewing motion. However, at 
high rotational speeds the roller end and the flange are often subjected to a oepletion 
in the lubricant supply due to centrifugal effects. In such cases, the minute amount of 
lubricant available at the roller-end-flange conjunction might well represent an example 
of steady-state starvation. Starvation effects in hydrodynamically lubricated contacts 
are important also if one wishes to calculate the rolling and sliding resistance and/or 
traction encountered in ball and roller bearings (!]. In another example, the effect of 
restricting the lubricant to a roller bearing is seen experimentally and theoretically 
to reduce the amount of cage and roller slip [2]. The theoretical analysis was accom- 
plished by changing the location of the boundary where the pressure begins to builduf. 
and noting the effect on the hydrodynamic forces. Combining this with relative velocity 
expressions and equilibrium equations enabled the determination of the amount of cage 
and roller slip. 

The location of the inlet and exit boundries as well as the respective boundary 
conditions to be applied has been one of the most controversial issues concerning star- 
vation of hydrodynamic contacts. The issue cf the effect of the lubricant supply on the 
inlet boundary condition and its consequences to incipient pressure buildup began to 
materialize as a result of earlier studies applied to rigid cylinders |3,4]. Laucer [5] 
and Tipei 16) asserted an upstream 1 mit of tt,? fluid film where the pressure begins to 
rise as governed by the Reynolds equation. This limit according to Lauder is determined 
by applying reverse flow boundary conditions (i.e., u •= 3u/8y * 0). Tipei locates 
the upstream limit as defined by the line of centers of two bounded vortices that are 
observed for pure rolling. Both cases have been criticized because their analyses lead 
to one position of pressure buildup regardless of the oil supply (7) . Dowsor i8] 

Floberg [9,10] and most dramatically, Wedeven, Evans, and Cameron [11] provide experi- 
mental evidence supporting the idea that the location of incipient pressure rise is 
determined by the oil supply. Further Wedeven, et al., using a Grubin type of EKD 
analysis, obtained very good correlation between experiment and the theory of starva- 
tion effects by choosing the start of the pressure buildup to occur at the meniscus 
boundary. Oteri (12), using stream function analysis for rolling rigid cylinders, 
showed that incipient pressure rise occurs at the meniscus boundary even in the presence 
of reverse-flow conditions. In view of this work, starvation effects in machine element 
applications can be predicted and relied on with a greater degree of confidence. 

One of the more important manifestations of lubricant starvation is the reduction 
in film thickness. This topic has received a good deal of attention in the literature 
111, 13-21). With the exception of |13, 19-21), these references are applicable only to 
elastodydrodynam.ic situations. Most of the work concerned with rigid contacts has oeer. 
devoted to line contact applications (13,20,21). Dalmaz and Godet (19) analyze the 
effect of the inlet on the film reduction factor for a sphere against a plate. However, 
to the authors' knowledge, an effort that parallels that of Hamrock and Dowson (14, lb] 
for the EilD contact is absent from the rigio contact theory. In those works, an 
expression was determined that relates the film reduction to the inlet cistance. 

The current study is a resumption of a previous r ig io-contact analysis 122] to 
extend validity for the minimum- f i Im- thickness equation derived there over a wider range 
of film thicknesses as well as to include the effects of starvation ir. this equation. 

The start of the pressure buildup as determined by the Reynolds equation is assumed to 
occur at the inlet meniscus. The location of the cavitation boundary was oetermineo ty 
applying the Reynolds boundary conditions as discussed in previous work [22]. The stuoy 
applies to a wide range of geometries (i.e., from a ball-on-plate configuration to a 
ball in a conforming groove). Seventy-four cases were used to numerically determine (1) 
an equation relating minimum film thickness with the fluid inlet level as well as with 
the dimensionless load-speed ratio and geometry, (2) an equation pr?dicting the onset of 
starvation, and (3) an equation predicting the onset of a critically starved con]unc- 
tion. The resulting equations are valid for dimensionless minimum film thicknesses 
Hp ranging from 5.0x10"^ to 1.0x10'^. Further, contour isobar plats and three- 
dimensional isometric pressure plots are presented. 

Numerical Procedure 


The hvdrodynamic effects on the central film thickness between two rigid solids in 
lubricated rolling and/or sliding contact are analyzed unoer conditions of lubricant 
starvation. The effects of starvation are determined by systematically cecreasing the 
fluid inlet level. The Reynolds boundary conditions are applied at the cavitation 
boundary, and zero pressure is stipulated at the meniscus or inlet boundary. The lubri- 
cant is assumed to be an incompressible Newtonian fluio under laminar, isotherm.al , 
isoviscous, and steady-state conditions. The numerical approach follows that of a 
previous investigation (22). There, a fully flooded film profile was specif lec and a 
pressure distribution satisfying the Reynolds equation was aetermined for a given speec, 
viscosity, and geometry. The analys’s treats the two rigid bodies as having parallel 
principal axes of inertia. This enables one to make a simplifying transformation to an 
equivalent system of a rigid solid near a plane separateo by a lubricant film (Fig. 1). 
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ReJcvant equations . - The same dimensionless expressions are used here as in |22]i 
that IE, 


a ] Eo 

The Reynolds equation 


X ■ x/R„, Y ■ y/Rv' ** ■ 


P “ pRx/ho“» ® • *y/Rx 


( 1 ) 


lx(“’ M) • Iv If) ■ I? 

is the governing equation within the conjunction. 


( 2 ) 


Kc recognize that, when the inlet supply levels art increased to values much 
greater than the minimum film thickness, calculations as governed by the Reynolds equa- 
tion are inherently in error far from the center of contact. The reason is that the 
Reynolds equation neglects curvature of the fluid film. Dowson |23J has pointed out 
that the errors involved in using this equation to determine the buildup of pressure 
in such regions are neglioible. The predicted pressures are themselves so very much 
smaller than the effective load-carrying pressures in the region of closest approach of 
the solids. The dimensionless film thickness equation is given as 


H + 1 
o 


- /I 


* a 


[l - - (Y/a)^] 


(3) 


where H is bounded above by the dimensionless fluid inlet level 


Ho (i.e. 


H 


i «ip< 


the dimensionless minimum film thickness 

The fluid inlet level Hjo is made dimeusAunx^ss ou knob, a,. ..h 
filled, Hjn * ^ conjunction is said to be fully flooded. 


Hin 

and below 

by 

e.g. 

, (Fig. 

1) . 


inle 

t is complet 

ely 

1 » 0 

at the 

cavi 

tat 

listr 

ibution 

that 



The Reynolds boundary conditions are used, that is, P * 3P/3N 
boundary and P » 0 at the inlet boundary (H = H^^) • R pressui 
satisfies these boundary conditions is then determined numerically by finite difler- 
encino with Gauss-Seidel successive over relaxation method. With this algorithm we are 
able to generate pressure distributions for given contact geometry, speeo, viscosity, 
film thickness, and fluid inlet level. 

Effect of inlet on pressure . - Figure 2 graphically portrays how a pressure dis- 
tribution is affected by the lubricant supply for a dimensionless minimum film thickness 
of 1x10"^. The two views, (a) and (b) , are for a fully flooded contact (Hm “ 1*00) 
and a starved contact (Hjn ' 0.001). The starvation effect on the pressure distribu- 
tion might be considered rather moderate. That is to say, the effect of reducing the 
fluid inlet level has essentially had no effect on the magnitude of the peak pressure. 
Only the size of the load bearing region has been affected. The load bearing region is 
defined by the meniscus and the cavitation boundary. 

We repeat this process of decreasing the fluid inlet level tor a given minimum film 
thickness for a wide range of minimum film thickness values. Figure 3 illustrates a 
situation in which the starvation is more severe. The minimum film thickness here is 
ten times as thick as what is indicated in Fig. 2. Note that the pressure distribution 
of Fig. 3 is not as localized as in Fig. 2. Consequently decreasing the size of the 
load bearing region has a more noticeable effect on the load capacity. Starvation is 
more evident here since we also set- that the peak pressure is noticeably reduced Iron 
the fully flooded condition. This sort of behavior takesplace if the fluid inlet level 
is of the same order of magnitude as the minimum film thickness. 

Generalized film thickness formula . - Now we've seen how starvation effects the 
pressure distribution for two different minimum film thickness values. In this analy- 
sis. we repeated this determination for minimum film thickness values ranging from. 

10'^ to 10"^. In all, 74 different cases were analyzed to determine a minimum film 
thickness equation as a function of dimensionless load to speed ratio (W/U) and fluid 
inlet level. 

The integration of the pressure distribution can be used in relating the hydro- 
dynamic effects (i.e., load, speed, and viscosity) to the minimum (central) film thick- 
ness for a given fluid inlet level. In general, for a given Hq snd u the load 
capacity w ano/cr the dimensionless load-speed ratio is determined as follows: 


n uR 
o X 


// 


P dx dY 


( 4 ) 


W/U 


// 


P dX dY 


where W/U • wnpURj^. 
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Kapitza 127), 
relationship: 


uein^ half Fommerfeld boundary conditions derived the following 


h 

V 


8»L VTu 


( 5 ) 


where L ■ n/2 half-Somerfeld B.C. (Ref. 27) 

L - 0.131 tan*^ (a/2) * 1.683 Reynolds B.C. (Ref. 22) 


In a previous work |22| , we found that if the problem is considered using Reynolds 
boundary conditions, then L required the above modification. But even then the equa- 
tion was valid only if the pressure distribution was very localized in the case of thin 
films (i.e. 5x10"^ 1 Ho ^ 10"'). Consider, for example, a ball in rolling motion 
that is loaded against a flat plate. For a load-speed ratio W/U of 340, the numeri- 
cally determined value (i.e., as determined by finite difference analysis) of the 
dimensionless film thickness is 10*^ (Table I). Equation (5) predicts the dimension- 
less film thickness to be 1.22x10'^, which is in error of the numerical value by 
22 percent. Consequently, we wish to revise Eq. (5) so that it is valid for the thicker 
films as well. This revised equation should reduce to Eq. (5) in the limit for thin 
films. Furthermore, the revised dimensionless film thickness should be expressed in 
such a way as to easily include the effects of starvation. This would then enable us to 
present one general exo't'.sion to be presented for tlie dimensionless film thickness that 
can be used for the f.J' range of film thicknesses for a starved conjunction as well as 
for a fully flooded c ” ,unct ii^n . Afttr the numerical analysis for each case was com- 
plete, the several curve fits or regression curves of Hq on W/U were considered 
that would be consistent with these above requirements. The most suitable curve fit 
considered was in the form of a more general linear equation, that is. 





( 6 ) 


Generalized film thickness formula (applicable to starveo as well as fully flooded 
conjunctions) . - The effect of lubricant starvation on the hydrodynamic film thickness 
was observed by varying the fluid inlet level to the contact ana noting the effect on 
load capacity for five different film thicknesses for the ball-on-plate contact (i.e., 
a « 1.00). In addition to the fully flooded data, 55 computer-generatea aata points 
were used to arrive at a family of equations having the form given in Eq. (6). An equa- 
tion for each fluid inlet level was determined by performing a linear regression by the 
method of least squares. Table 11 lists for each fluid inlet level the values of the 
coefficients Cj and Cq, the coefficient of determination r’, and D;nax' maximur. 

percentage of error D defined as 


D 



X 100 


( 7 ) 


Note that Cj remains essentially unchanged, and can be determined from 


1 

4L(128a) 


T7I 


18 ) 


iRat 


Furthermore, all of the effect of starvation is described in the value of Ar 

expression for the coefficient Cq as a function of the fluid inlet level H 
would enable a determination of a generalized minimum- f i Im-thickness formula 
applies to starved as well as fully flooded conditions. A close eramination of the 
variation of Cq with Hm in Table II reveals that Cq« ( 2/H jp) for the severely 
starved situations. As approaches 1, Cq approaches a value very neatly equal 

to e (i.e., the base of the natural systems of logarithms, 2.718). This suggests 
e^in as a modulating factor. Further considerations for the nearly floooed inlet 
levels show that 
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Finally for the full range of values, Cq varies with Hm as follows: 

1/2 


1.11 e 




(V) 


Thus, our generalized minimum film thickness formula in terms of geoit.etry (i.e., 
radius ratio a), load-speed ratio W/U, ana fluid inlet level Hjq, can be written as 


W/U 
4L (128a) 


T77 


♦ 1 




-2 


in 


(lOj 


7 he n'c.Tuir of aarpement between the calculated and input values of Hq ia re- 
pre^ertro hv the value of D (Eq. (7)) and presented in Table 1. Table I ahows that, 
vhrr the fluid inlet level is of the same order of magnitude as the minimum film thick* 
nesf , the error that rerults from using Eq. (10) becomes larger. However, Eq. (10) can 
confinertly be used for the full range of minimum film thickness if the fluid inlet 
level IF Fiich that 0.00< <, * 1.000. For very thin films (i.e., Hq * 10'^) Eq. (10) 

can t>e useful throughout the full range of fluid inlet levels that were~invest igated 
(i.e.. 0.001 < Hjp < 1.000). Note that the film thickness formula is intended to be 
used only for~a range of speeds and loads in which piezoviscous and deformation effects 
are negligible. It was determined that these effects can be significant for minimum 
film thicknesser less than 5.0x10"^. 


Note also from Table I that excellent agreement is obtained for the near*line* 
cootact applications (i.e., a ■ 36.54) even though these data were not used in the 
determination of the above Eq. (10). Host of che predictions by Eq. (10) are within 
3 percent of the numerically determined values and do not exceed 6 percent for any case. 


Fquatiop (10) is equally valid for the fully flooded conjunction as well as for the 
Etarver conjunction. That is, if we set Hm ■ 1.00, then Eq. (10) reduces to 



W/U 


vTTBa 


♦ 3.02 


( 11 ) 


fquatior (11) should be used in place of Eq. (5) whenever a calculation is made for 
a fully flooded conjunction since it is valio fer a broader range of film thicknesses 
.'j.e., to 1x10"^). Table 1 shows that the error for the full range of mini- 

rr.;iT filn thicknesses investigated (for Hjf, ■ 1) is less that 1 percent for a ■ 1. 

Ff r the near-line-contact geometry (i.e., a « 36.54) the error does not exceeo 3.36 
percent . 


Reduction in minimum film thickness . - It is now possible to determine the reduc- 
tmr in minimum film thickness from the fully flooded value if the fluid inlet level is 
known. This can be done by inserting Eq. (11) into Eq. (10). 



Civicir.g both siner of the equation by 


Hq gives 



( 12 ) 
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(13) 


wricrn t IF the reduction in minimum film thickness due to starvation. 


PFSL’I':.': ANr DISCUSSION 


Ff'ert of starvation 

on 

pressure distribution. 

- The discussion of lubricant star 

vatjor car be facilitated 


focusing on one of the 

simplest geometric arrangements 


(I.e., a ball rolling and/or sliding against a flat plate) as shown in Fig. 2. The 1 ig- 
ure comparer the pressure distribution determined numerically for the f ul ly- 1 loooec; 
inlft with the most severely starved inlet (Hm “ O.OCl). The comparison is mace for 
a constant minimum^ film thickness (i.e., Hq “ 1.0x10"^). Note that the pressure 
peak built up in the starved inlet is only slightly smaller than that ot the fully 
flooded inlet. However, the area of pressure build-up is considerably smaller, and sc 
the starved inlet is unable to support as much load for a given film thickness as the 
f ul 1 y- f 1 ooded inlet. 

Fiacre T provioes the same sort of comparison but for a thicker miinimum film (i.e., 
Hf) * 1.0''10*'). The significant difference between the two figures is that the 
starved inlet for the thicker film has a more pronounced effect or the pressure peak. 

The fluid inlet level (Hjn ■ 0.002) for the thicker film represents a relatively more 
hiahlv starved inlet since Hjp is of the order of in this case. The other 

featu e to be noticed in comparing Figs. 2 and 3 is that the pressure oistribution is 

more evenly spread out for the thicker film. Thus, changes in the meniscus (or integra- 
tion domain) are going to have a more noticeable effect on the load-carrying capacity. 
Note also that because of the boundary conditions the integration domain takes on a 
"kidney-shaped* appearance. This is more clearly shown in the isobaric contour plot 
shown in Fig. 4. 

Mini m um film thickness er^uation . - Thus far we have compared the pressure build up 

in a sevefely starved inlet with that in a fully floodeo inlet for two minimum film 

thicknesses. Our investigation, however, included several fluid inlet levels for a 
variety of minimum film thicknesses. The results are summ.arized in Table 1. 
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A 9 eneralited minimum film rhicknets formula (Eq. (10)) wai derived from the 
reaults of Table 1. Fiqure 5 indicates how well the equation represents the computer 
generated data in the table. It was not possible to display all these results in Fig. 

5. However, the figure is representative of the overall results. The equation fits the 
data quite well excep*^ when the fluid inlet level is of the same oroer of magnitude >s 
the minimum film thickness. Based on the discussion concerning peak pressure, it woulo 
seem that ‘.he formula holds well for those cases in which the degree o' starvation is 
such that peak pressure is not significantly reduced. 


Of course, it would be most desirable to compare the data in Table I with experi- 
mental data. To the authors' knowledge, the only available experimental data were 
obtained by Dalmaz and Godet [25]. To compare Eq. (10) with experiment, the data from 
125] were replctted in Fig. 6. The experimental data were taken under lightly loaded 
(rigid contact), isoviscous conditions for pure sliding of a ball on a plate. The fluid 
inlet level in these experiments was reported to be 1 millimeter. The '^all diameter was 
30 millimeters; consequently the dimensionless fluid inlet level 0.067. The 

experimental data were presented as a plot of the dimensionless parameters Hq/WG 

versus The materials parameter G in the plot was included so as to ODIAIMAl DBAr ' ia 

accommodate the elastohydrodynamic range in a more general way. rAvit fS 

OF POOR OUALfTr 

Here, we wish to compare our hydrodynamic starvation theory only with the hydro- 
dynamic results of Dalmaz and Godet. To do this, the ordinate ana abscissa were proper- 
ly scaled (assuming a reasonable value of WG ■ 4.5x10*^) so that the minimum film 
thickness could be plotted against the dimensionless load-speea ratio as shown in 
Fio. 6. The solid line in Fig. 6 is a plot of Eq. (10) for a » 1 and « 0.067. 

The dashed line represents a previous theory |22] in which the reouced inlet level due 
to starvation was not considered. The present theory (Eq. (10)) shows better agreement 
with the experimental results than our previous theory |22] for lower values of K/D. 

This is because under conditions of low load and/or high speed, the pressures become 
more significant away from the center of tne contact. Consequently, neglecting the size 
of the inlet introduces an increasing amount of error as K/u is decreased. 


Although consideration of the size of the inlet domain improves the agreement of 
the theory with experiment for the thicker films, still further improvement is possi- 
ble. It is believed that the effects of reverse flow in the inlet must be includeo in 
the theory to obtain better agreement. If reverse flow is considered, not all the 
available lubricant determined by the fluid inlet level will pass through the contact. 
The hydrodynamic contact would essentially s?e this as a reduction in supply from what 
actually is there. In other words, the inlet is more severely starved than we have 
taken into account. From Fig. 5, we see that increasing the severity of starvation has 
more effect on the load capacity for thicker films (i.e., Hq « 10*-') than it ooes 
for thinner films (i.e.. Ho “ 10"^). Consequently, reverse-flow considerations 
should improve agreement between theory and experiment for the thicker films while still 
maintaining good agreement in the thin-film range. 


Lubricant film thickners reduction factor and onset of starvation . - Of practical 
importance to luf^'icant starvation is the reduction in minimum film thickness from the 
fully flooded value. Equation (13) is a derived expression for & in terms of the 
fluid inlet level and the fully flooded film thickness (also given by Eq. (11)). 

Figure 7 is a plot of B as a function of the fluid inlet level H^p for several 
values of f. It is of interest to determine a fully f loooea- star veo bouroary 

(i.e., that fluid inlet level after which any further decrease causes a significant 
reduction in the film thickness). Hamrock and Dowson |14] oetermineo this bouncoiy for 
elastohydrodynamic (EHD) applications upon satisfying the following ccnoitior.: 


1 


6 , « 0.03 

”in • 


(14) 


The value of 0.03 was used in Eq. (14) since it was ascertaineo that the oata ir. Table I 
were accurate to only ±3 percent. 


Thus, for a given value of )^p,f* can solve for a value of 

satisfies Eq. (14). A suitable relationship between H*p and 
tamed by generating a table of values (e.g.. Table III) and fitting 
the method of least squares. This gives 


hJp that 
can be ot>- 
a power curve 


by 


H,„ - 4.11 (H^ .) 
in o , r 


0.36 


(15) 




i 

> 


Thus, w^ have an equation that determines the onset of starvation. That is, for 
Hin i Hin • fully flooded condition exists, whereas for Hjp < Hm a starvec 
condition exists. 

Critically starved inlet . - In certain bearing applications, the power loss result 
ing from churning of the oil may be higher than the power loss resulting from friction 
of the bearing alone 126). These power losses can be minimized by reducing the lubri- 
cant supply until a loss in film thickness causes the friction losses to increase. 
According to the results shown in Fig. 7, the fluid inlet level car. be decreased sub- 
stantially without adversely affecting tlie minimum film thickness. Consequently, it 
might be advantageous to operate the bearing with a lubricant supply ]ust sufficient to 
preclude any drastic reductions in minimum film thickness, as seen in the figure. Such 
a critical fluid inlet level might well be defined to occur at the knee of the curve, 
that is. 
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- 1 (Ifc) 

^in “ *^in,cr 

Solutions to this expression for several values of ininimuir film thickness are listeo in 
Table III. A power curve fit by the method of least squares gives 

0.2977 

"m.cr - 0.875 (H^^^)“*^’” (17) 

The regression coefficient for this fit was aetermined to be 0.9999, inoicating an 
extremely good fit. 

CONCLUDING REMARKS 

Numerical methods were used to determine the effects of lubricant starvation on the 
minimum film thickness under conditions of hydrodynamic point contact. Starvation was 
effected by varying the fluid inlet level. The Reynolds bounoaiy conoitions were 
applied at the caviation boundary, and zero pressure was stipulateo at the meniscus or 
inlet boundary. The analysis is considered valid for a range of speeds and 1 ~>ads for 
which thermal, piezoviscous , and deformation effects are negligiole. It can he applieo 
tc a wide range of geometries (i.e., from a ball-on-plate configuration to a ball in a 
conforming groove). Seventy-four cases were used to numerically determine 

(1) A generalized expression for the minimum film thickness as a function of cimer- 
sionlesE load-speed ratio, geometry, and fluid inlet level (Lq. (10)). The expressior 
should be applied for film thicknesses in the range 1.0x10*^ i f^o i l.OxlO'^ ano for 
fluid inlet levels of 0.004 < Hjn ^ 1.00. For 5x10"^ i ^>o 1 10"^ the equation 
can be applied for a fluid inlet level of 0.001 < Hjp ^ 1.00. 

(?) A film thickness reduction factor (Eq. (13)) expressed as a function of tf.f 
degree of starvation (or fl^iid inlet level) for a given fully flooded film thicKness 
vaiue . 

(3) An equation (Eq. (IS)) that determines the onset of starvation. 

(4) An equation (Eq. (17)) that determines a critically starver contact. Contour 
I'obar plots and three-dimensional isometric plots also pi ‘senteci. 
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TABLE I*. - COWARISOB OF DIHCNSIOHLESS FILN THICOIESS mUES FROM EQUATIOM (10) 
WITH OIWRSIORLESS FILM THICKMESS VALUES INPUT TO COMPUTER 
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TABLE m*. - OITCNSIUNUSS FLUID INLET VALUL:> THAT DITIRHINE STARVED - 
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VALUES OF HINIHUN FILM THICKNESS FOR A FLOODED LUNJUNCTIUN 
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Figure 1. - Depiction of lubricant flotw for a roiling/sliding contact and 
corresponding pressure build-up. 


Rolling 



I 



(b) Starved condition: dimensionless fluid inlet level Hj^-UGOl. 

Figure Z - Three-dimensional refresentation of pressure distribution, comparing starved with fully 
flooded conjunction for dimensionless minimum film thickness • 1x10"^. 
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(a) Fully flooded condition: 
dimensionless fluid inlet level 
Hjp - 1.00; dimensionless 
maximum pressure • 

1. 20Kl(r; dimensionless load- 
speed ratio W/U- 1133L6. 


(b) Starved condition: 
Hin-aWP™x-1.19xl(^: 

m-86Z6. 


;c) starved condition: 
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Figure 4 - Isobaric contour plots for three fluid inlet levels for dimensionless minimum film thickness 
Hg - IxlO"^. 
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o 1.000 
O .035 

o .010 

□ .004 

A .002 

o .001 

Eq. (10) 
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Figure 5. - Comparison of dimensionless minimum film thickness equation (eq. (10)) 
with computer-generated data as a function of dimensionless load-speed ratio 
several values of dimensionless fluid inlet level. 
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Figure 6. - Theoretical and experimental results. 
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Figure 7. - Minimum film thickness reduction factor as a , 

function of fiuid inlet level for several values of dimen- 
sionless minimum fiim thickness for flooded conjunction. 
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